Dynamics of a tight-binding ring threaded by time-periodic magnetic flux 



W.H. Hu and Z. Son£| 

School of Physics, Nankai University, Tianjin 300071, 



China 



o 

(N 

> 

O 



43 : 

Oh- 
i , 

> , 

c : 
a . 
3 ■ 
a* 



> 

in 
o 



We analytically study the effects of periodically alternating magnetic fields on the dynamics of a 
tight-binding ring. It is shown that an arbitrary quantum state can be frozen coherently at will by 
the very frequent square- wave field as well as the monochromatic- wave field when the corresponding 
optimal amplitudes are taken. Numerical simulations show that the average fidelity depends on 
not only the system parameters, but also the features of the quantum state. Moreover, taking the 
initial zero-momentum Gaussian wave packets as examples, we show the dependence of the threshold 
frequency on the width of the packet for the given average fidelities. These observations provide a 
means to perform the quantum state engineering. 



PACS numbers: 03.67.-a, 03.75.-b, 03.75.Lm 



I. INTRODUCTION 

Coherent quantum state storage and transfer via a 
coupled qubit system are an important problem in the 
emerging area of quantum information processing (QIP). 
One of the promising method of quantum state transfer 
is employing a solid-state data bus with minimal spa- 
tial and dynamical control over the on-chip interactions 
between qubits [3-01- However, small imperfections in 
receiving a quantum state and storing it coherently can 
seriously affect the fidelity of QIP. Stopping and freez- 
ing a flying qubit within a region of the data bus is a 
tool for this task. It has been proposed that a cou- 
pled cavity array system exhibits the possibility of an 
all-optical coherent control of light [81-KLlj. The dynam- 
ical control includes adiabatic scheme, under which the 
quantum state is fixed on the superposition of the instan- 
taneous eigenstates when the Hamiltonian varies slowly, 
and bang-bang control techniques, by means of a dynam- 
ical control field, at averaging to cease the unwanted evo- 
lution of the state. 

Dynamical decoupling (DD) is a well established 
paradigm of bang-bang control techniques, which em- 
ploys a specially designed sequence of control pulses ap- 
plied to the qubits in order to negate the coupling of the 
central spins to their environment 



12|. Moreover, the 



quantum Zeno effect has been proposed as a strategy to 
protect coherence [H, HH ■ Recently, it is proposed that 
a periodically driven potential can suppress the tunnel- 
ing between adjacent sites in a lattice [TH, [l6j]. In this 
paper we will pay attention to a fundamental aspect of 
QIP and generally study the influence of periodically al- 
ternating magnetic fields on the dynamics of a quantum 
state on a tight-binding ring. We consider the dynamics 
of the states in a tight-binding ring system that is pierced 
by a time-periodic magnetic flux $ it) = $o + $a/ (wt) 
with the angular frequency u>. We investigate the impact 
that the amplitude and frequency of the flux might have 
on the efficacy of the quantum control. We focus on the 
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suppression of the evolution through a bang-bang control 
procedure, and study how the occurrence of a controlling 
field modifies the effectiveness of the control procedure. 
Our analysis is focused on the behavior of the fidelity of 
the evolved state with respect to the initial state, which 
has been employed to measure the efficiency of quantum 
state transfer. We will show that the time evolution of 
a state in such a time-dependent Hamiltonian can then 
be treated as an adiabatic process without any approxi- 
mation. Then analytical results can be obtained, which 
should give more insight into quantum measurement and 
control. Additionally, this scheme can be applied to a 
neutral-particle system by torsional oscillation. 

In Sec. |TI] we derive a general formalism for such a 
time-dependent system. We further introduce the ex- 
pressions for the time-averaged fidelity, thus completing 
the description of the controllability of the systems for 
the quantum states. Sec- IIIIllTVl are devoted to the appli- 
cations of the formalism. These include a detailed treat- 
ment and computation of the time evolutions of typical 
initial states under the square and monochromatic time- 
periodic flux, respectively. Final conclusions and discus- 
sions are drawn in Sec. [Vj 



II. MODEL AND GENERAL FORMALISM 

In this section, we present the charged particle model 
under consideration, a simple tight-binding model in an 
external magnetic field. Here, the particle-particle inter- 
action is ignored for simplicity. Our approach is based 
on our previous work in Ref. [17| . where we have pro- 
posed a scheme for quantum state transfer. It employed 
a loop enclosing a static magnetic flux to control the 
speed of a wave packet. Another basic operation for the 
quantum-state engineering is coherently freezing a state 
on demand. For instance, quantum information process- 
ing requires transferring, stopping and freezing a flying 
qubit within a region of the data bus. In Ref. [l7| . we 
studied how to move a Gaussian wave packet at a certain 
speed on demand by a static magnetic flux. In this work 
we aim at employing the same system with a periodi- 
cally alternating flux for freezing a wave packet. Com- 



2 




• — •- 



<D(t) 



• » • -» 




where k = 2im/N, n S [1 , iV] , the Hamiltonian can be 
readily written as 



H = -2 cos ( k + <t> (*)) a W, 



(3) 



ft 



with </> (i) = 27r<E> (t) /AT and the corresponding eigen- 
states in the form of 



li>. 



(4) 




FIG. 1. (Color online) The schematic illustration for a tight- 
binding ring threaded by a time-periodic magnetic field. The 
time-dependent flux for demonstrating the influence on the 
evolution of a quantum state is in two types: a square and 
a monochromatic waves, respectively, with amplitude (j>A and 
period r. 



bining the two schemes, one can accomplish the task of 
"coherently storage-transfer-storage" of a quantum state. 
Now we will generalize this description of the system in 
Ref. fl7| by allowing for an additional time-periodic flux. 
We restrict our attention to the influence of the applied 
periodically alternating field on the dynamics of the par- 
ticles. 

Consider a ring lattice with N sites threaded by a mag- 
netic field illustrated schematically in Fig. [T] The Hamil- 
tonian of the corresponding tight-binding model 
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H(t) = -Jj2( 
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e «***(t)/N a i a . +1 
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(1) 



depends on the magnetic flux through the ring in units 
of the flux quantum $o = hje. Here at is the creation 
operator of a particle at the jth site with the periodic 
boundary condition cin + i = a\. The flux does not exert 
force on the Bloch electron, but can change the local 
phase of its wave function due to the Aharonov-Bohm 
(AB) effect. Note that the particle is not restricted to be 
cither fcrmion or boson. 
Taking the transformation 
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e ikj a k7 
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Note that the time-dependent Hamiltonian possesses 
fixed eigenstates, while the flux solely affects the eigen- 
values. This will be crucial to employ such kind of setup 
to investigate the control of a quantum state due to 
the rareness of the exact solutions to a time-dependent 
Hamiltonian. The evolution of an arbitrary state under 
the Hamiltonian H is dictated by the unitary operator 

U (i', t) = exp C-if?Hdt"\ which yields the propaga- 
tor represented in the momentum and spatial eigenstates 
as 



U k , k {t',t) = (k'\U(t',t)\k) 
U fj (t',t) = {j'\U(t',t) \j) 

iV f 't(j'-j) e i2Jft(''f) 



(5) 



where 



fk(t',t) 



cos [k + (j> (*")] dt'' 



(6) 



We note that the propagator is in diagonal form in k 
space. 

In general, one employs the fidelity 



F(t) = |ty(0)|CT(t,0)hKO)>| 



(7) 



to characterize the relation between the target state and 
the evolved state at time t. However, when ip (t) is a 
periodic function, F (t) should be oscillating. Thus the 
long time average will be appropriate as a measure for 
the deviation from the original state. 
The average fidelity is defined as 



F 
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lim — 



F (t) dt 



(8) 



dt, 



where 



c k = (k |V (0)) 



In this paper, we focus on the case of periodic (j> (t) with 
a period of r 



0(t) = 0(t + r). 



(9) 
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In the following sections, we will apply the general for- 
malism to the cases of square and monochromatic waves 
due to the following reasons: The square- wave case is an 
demonstrative example since it is the simplest model to 
calculate. The exact solution for this particular case is 
helpful to clearly present the main idea of the scheme 
without involving much more complicated calculation. 
And the monochromatic-wave case is a more practical 
situation. We will argue that the quantum state freezing 
can be achieved by the high frequency alternating flux 
when its amplitudes $a are optimal. 



III. SQUARE WAVE 

Let us begin the discussion with the simplest case: the 
flux is in the form of 



which leads to 



4> (t) = 4> + ^ASgn (sin (wt)) 



(10) 



where 10 = 2tt/t is the angular frequency and sgn indi- 
cates the sign function. It is a toy model which demon- 
strates how, in principle, the periodic alternating flux can 
prevent a quantum state from spreading. The Hamilto- 
nian can be diagonalized as 

H = -2J^ £fe (*)4a fe , (11) 

k 

£k {t) = cos [k + </>ASgn (sin (u)t))] (12) 

where we absorbed (f>o into k by k — > k — ipo for the sake 
of simplicity. Then we have 



k 



\ck\ 2 exp ( i2J / e k (t')dt' ) dt. 



F (t) = lim - 

(13) 

Obviously, it is hard to get the analytical expression of 
F (r). However, one can get insight into the influence of 
the square wave on a quantum state from the following 
analysis. During each interval, with the flux being static, 
the dynamics of the quantum states in such situation has 
been discussed in Ref. [l7j ]. In particular, for 0a = vr/2 , 
the corresponding Hamiltonians are 

H{t) = {u +1 Sgn S Sin H!!^n (14) 

w \ sgn (sin (wt)) -< v ; 

with 

H± = ±2J^smka\a k . (15) 

k 

Then the dynamics on the successive time intervals 
[t, t + t/2] and [t + t/2, t + t] are time reversal processes 
to each other, i.e., 

U k 'k (t + r/2, t) = U k , k (t + r, t) U k , k (t + r/2, t + t) 

= Uil{t + T,t + r/2), (16) 
*e [<W2], 



foHt)> = hK« + r)) 



(17) 



for an arbitrary state. Then after a period of time r, any 
state will go back to its initial state, i.e., 



F{t)=F(t + r), F(nr) = 1. 



(18) 



For small t, the evolved state should not leave its initial 
state so far at any time. Intuitively, in the limit of r — > 
0, the initial state may be frozen at its initial position. 
Actually, the periodicity of \ip (t)) admits 
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Y,\ck\ 2 e l2JsiDkt 



dt, 



(19) 



and 



lim F ~ lim — 



r/2 



1 + i2Jt \ c k \ 2 sink 



dt = 1, 



due to J2 k \ c k\ = 1, 



(20) 

^ fc |cfe| 2 sinfc -< 1. We conclude 

that an arbitrary state can be frozen at will when the 
frequency is sufficient high. Particularly, in the case of 
\c k \ being symmetrical about k = 0, we simply have 

r/2 



F= 2 - 

T 







|cfc| 2 cos (2Jt sinfe) 



k 



dt. 



and for small r (r <C J 1 
1 



it has the form of 



F^l-ijV£| Cfc | 2 sm 2 fc. 



(21) 



(22) 



Obviously, ^ fc |cfc| 2 sin 2 k determines the speed of con- 
vergence: the more the contribution of \c k \ 2 to zero k, 
the faster the convergence speed. In the following, we 
estimate the convergence speed of a specific state. A 
Gaussian wave packet with central momentum ko can be 
expressed in the form of 



Cfe = A exp 



ko) 



(23) 



where 



A = 



exp (—a 2 (k — fco) 2 ) 



is the normalization factor and equal to yj N / (2y / 7ra 2 ) 
for large N. Here a determines the width of the wave 
packet. We consider such a state with fc = as the 
initial state. For large frequency v — 1/r, we have the 
frequency dependence of the average fidelity as 



F — 1 



J 



12^ 2 



(24) 
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FIG. 2. (Color online) Average fidelity as a function of field 
amplitude 4>a (units of it) and frequency v (units of J) of a 
square-wave flux for the Gaussian wave packet with a = 50 
and central momenta ko = 0, n/4, i/2 on an JV = 1000 
ring. The average fidelity is computed over the time interval 
[0, 25N/J]. It shows that tt/2 is the optimal amplitude and 
the fidelity becomes very sensitive to the magnitude for the 
initial wave packet with high speed. 



FIG. 3. (Color online) Average fidelity as a function of fre- 
quency of the square-wave flux for a Gaussian wave packet 
with q = 50 on an N — 1000 ring. Panels stand for different 
values of the amplitudes 4>a = tt/8, 7r/4, 3tt/8, tt/2, 5n/8, 
37r/4, 77t/8 and central momenta ko = 0, 7r/4, 7r/2 of the ini- 
tial wave packet. It shows that, when 4>a = tt/2, frequently al- 
ternating flux suppresses the evolution of the quantum states. 
Here the frequency v is expressed in units of J. 



Then for the given initial state, we have evaluated the 
threshold frequency as 

v c = J J 1 T e ~° , (25) 
V 12(1-3^' v ' 

from which one can obtain the control of the quan- 
tum state with fidelity F c (r) . We plot Eq. (1231) in 
Fig. 03 which shows that the threshold frequency in- 
creases rapidly when the width of the state gets narrow. 

So far, our investigation is carried out analytically for 
a particular case of </>a = 7r/2. It can be seen that 
the mechanism of the perfect quantum state freezing ex- 
pressed as Eq. (|T8|) is the periodicity of the evolved wave 
function Eq. (|17j) arising from the amplitude 4>a = tt/2. 
It is worthy to investigate what happens for other values 
of 4>a and finite frequencies. In the following, we will 
perform a series of numerical simulations from several 
aspects. 

First of all, we investigate the influence of the ampli- 



tude. The numerical simulations are performed in Fig. [2] 
for a Gaussian wave packet which has the form of Eq. |23|) 
with different central momenta in the systems with dif- 
ferent amplitudes and different frequencies. It shows that 
the average fidelity approaches to unit for (/>a = 7r/2 
when v is sufficiently high, which is in agreement with 
the above analysis. And it is noted that the locations of 
the maxima of the average fidelity for moderate v shift 
to the left for the case of ko = 0. It is important for 
practical implementation. In the case of not sufficiently 
high frequency, the optimal amplitude should be smaller 
than 7r/2. For ko = 7r/4, 7t/2 the fidelity becomes very 
sensitive to the magnitude of 4>a- This feature can be 
exploited to select the wave packet with preferable fco in 
the following way. Suppose a many-particle initial state, 
which consists of wave packets with various speeds. One 
can first tune 4>o t° meet 4>q + ko = tt/2. Then take 
4>x = tt/2 to hold the wave packets with ko on demand. 

Secondly, the analysis above is based on the assump- 
tion that not only the amplitude of the flux is 7r/2, but 
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FIG. 4. (Color online) Average fidelity as a function of fre- 
quency of the square-wave flux for a single-site state on an 
N = 1000 ring. Plots are presented for different values of the 
amplitudes 4>a = tt/8, n/4, 3n/8, n/2, 5n/8, 3ty/4 and 7n/8. 
It shows that, when (j>A = t/2, frequently alternating flux 
suppresses the evolution of the quantum states, while other 
cases exhibit small fidelity. Here the frequency v is expressed 
in units of J. 



also the frequency is sufficient high. Then even the ampli- 
tude is taken exactly as 7r/2, the efficiency of the scheme 
is different for different quantum states under the finite 
frequency. For instance, the average fidelity during the 
period r mainly depends on the overlap of the initial 
wave packet and its evolution driven by H±. In order 
to demonstrate these analysis, the numerical simulations 
are performed for two kinds of initial states: a Gaussian 
wave packet with central momentum kg, which has the 
form of Eq. (|23l) . and single-site state |Z) = at |0). We 
consider the time evolutions of this GWP with different 
fco in the system with different amplitudes and different 
frequencies. The average fidelity F (r) over the interval 
T ^ 25N/J is plotted in Fig. H It shows the follow- 
ing features: (i) The average fidelity approaches to unit 
for all the given initial wave packets with different fco 
when the amplitude is 7r/2, which is in agreement with 
the above analysis, (ii) The threshold frequency in the 
case of 4>a — tt/2 gets lower as fco goes closer to 7r/2. 
This also accords with the above analysis. Actually, the 
velocity of a 7r/2 wave packet becomes zero under the 
Hamiltonians H± . Thus it deviates from the initial state 
slightly during the interval t, leading to a high fidelity, 
(iii) The optimal average fidelity becomes more sensitive 
as fco goes closer to 7r/2, which is in agreement with the 
results in Fig. [3J 

Based on these features, one can design a scheme to 
achieve the maximal fidelity. In ideal case, for any given 
frequency of the field, </>a = tt/2 is always preferable. 
When the frequency is not sufficient high, one can tune 
cj)Q to match fco in order to achieve a lower threshold fre- 
quency. However, the accuracy of the field may affect the 
fidelity due to the sensitivity of it around 4>a — tt/2 in 
practice. Then one can tune 4>q to stabilize the fidelity. 



FIG. 5. (Color online) Threshold frequency on the widths 
of the zero- momentum Gaussian wave packet. The plots are 
Eq. Q25p and numerical results corresponding to the average 
fidelities as 0.96 ± 0.01, 0.90 ± 0.01 and 0.85 ± 0.01. It shows 
that the threshold frequency increases rapidly when the width 
of the state gets narrow. 



Finally, in order to demonstrate the applicability of 
our findings to control a quantum state, we also plot the 
average fidelity for a single-site state \l) = a, |0) in Fig. 2) 
Such a state has \ck\ 2 = l/N and is the narrowest limit 
of a wave packet. It shows that the average fidelity ap- 
proaches to unit only in the case of <pA = t/2 and a rel- 
ative high frequency. And to demonstrate the efficiency 
of this method, we take a initial zero-momentum Gaus- 
sian wave packet as examples. Numerical simulation is 
performed in Fig. [5] for the dependence of the threshold 
frequency on the widths for the given average fidelities. 
For a comparison we draw the curves from Eq. (|25[) and 
numerical results corresponding to the average fidelities 
as 0.96±0.01, 0.90±0.01 and 0.85±0.01. The width coef- 
ficient a of the wave packet is taken from 1 to 50, ranging 
from a single-site state to a very wide wave packet which 
is approximately a plane wave with zero momentum. As 
can be seen from the figure, analytical and numerical re- 
sults both indicate that the threshold frequency increases 
rapidly when the width of the state gets narrow. Then 
for a finite-frequency field, such a scheme has a high ef- 
ficiency for a broad wave packet. 



IV. MONOCHROMATIC WAVE 

The simplicity of the square- wave field makes it easy to 
take an analytical investigation for the concerned prob- 
lems, because the exact solution of this model is helpful 
to clearly present the main idea of the scheme without 
involving much more complicated calculation. However, 
such a toy model is not exactly accessible in experiments 
due to the sudden change of the flux. In this section, 
we will consider the monochromatic-wave field, which is 
more practical. It will been shown analytically and nu- 
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FIG. 6. (Color online) Average fidelity as a function of field 
amplitude </>a (units of n) and frequency v (units of J) of 
a monochromatic-wave flux for the Gaussian wave packet 
with a = 50 and central momenta ko = 0, n/4, tt/2, on 
an N = 1000 ring. The average fidelity is computed over the 
time interval [0, 25N/ J] . It shows that 0.7657T is the opti- 
mal amplitude and the fidelity becomes very sensitive to the 
magnitude for the initial wave packet with high speed. 



The monochromatic-wave field is in the form of 
4> (t) = 4>q + 4>a sin ut. 



(26) 



Unlike the square-wave, even in the special case of </>a = 
7r/2, the fidelity is no longer a periodic function due to the 
breaking of the time reversal symmetry. Thus one should 
consider the integral to the whole time duration. How- 
ever, the analytical function of <fi (t) may lead to some 
analytical results. Here we still neglect </>o for simplic- 
ity. We firstly investigate some special cases analytically 
to seek the optimal 0a satisfying the relation Eq. ([T5]). 
and then perform numerical simulations for more general 
cases. 

Considering the evolved state at the instant nr, where 
n is an integer, we have 

r-nr 

fk(nr)= / cos [k + <px sin lot] dt' (27) 
Jo 

— nr cos k Jo (4>a) ■ 



FIG. 7. (Color online) Average fidelity as a function of fre- 
quency of the monochromatic-wave flux for a Gaussian wave 
packet with a = 50 on an N — 1000 ring. Panels stand for 
different values of the amplitudes <f>A = 7r/8, 7r/4, 37r/8, n/2, 
57r/8, 0.765-7T, 77r/8 and central momenta ko = 0, 7r/4, tt/2 of 
the initial wave packet. It shows that, when <f>A = 0.765-7T, fre- 
quently alternating flux suppresses the evolution of the quan- 
tum states. Here the frequency v is expressed in units of J. 



Here 

i r 

J m (aO = - / cos {m6 ~xsin(9))d9, (28) 
t Jo 

are the Bessel functions of the first kind. Then the cor- 
responding fidelity is 



F(nr) 



E 

k 



| C J 2 gi2.7nrcos kJ a (<f> A ) 



(29) 



Note that taking <pA = 0.765-7T, we have j7o (4>a) = 0, 
which leads to F (nr) = 1 for an arbitrary initial state. 
This fact is quite similar to the case of the square-wave 
field with cf>A = tt/2. 

Now we turn on our numerical investigation to the 
more general cases. We perform the numerical simula- 
tions for the same states discussed in the last section. The 
numerical results are plotted in Figs. |H [7] and [5J They 
show that the square and monochromatic waves lead to 
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FIG. 8. (Color online) Average fidelity as a function of fre- 
quency of the monochromatic-wave flux for a single-site state 
on an TV = 1000 ring. Plots are presented for different values 
of the amplitudes (j> A = n/8, tt/4, 3tt/8, tt/2, 5tt/8, 0.765tt 
and 7n/8. It shows that, when cf>A = 0.7657T, frequently alter- 
nating flux suppresses the evolution of the quantum states, 
while other cases exhibit small fidelity. Here the frequency v 
is expressed in units of J. 
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FIG. 9. (Color online) The same as Fig. [5] but only numerical 
result for the monochromatic-wave field. 



above sections, we conclude that the evolution of a quan- 
tum state can be suppressed through the time-periodic 
flux. In both situations, the efficiency of the schemes 
depends on the parameters <pA, 4>q, ko and a in the 
similar manner. The features can be exploited to control 
quantum dynamics for quantum information and compu- 
tation purposes. 

V. DISCUSSION 

We have seen that the threaded magnetic flux, in stead 
of the electric field, plays an important role in controlling 
a state. It can be applied to a more extended system to 
control an uncharged particle. Actually, if the system is 
rotated, an effective magnetic field will be induced in the 
rotating frame of references. Therefore, a neutral parti- 
cle state in a ring lattice can be controlled via torsional 
oscillation. 

For a rotating ring with angular frequency f2, an addi- 
tional term 



-9.L, 



-QK 



N 



ia'a J+ i + H.c. 



(30) 



3=1 



should be added on the Hamiltonian with cf) = in the 
non-inertial frame [l8j], where if is a constant depends 
on the geometry of the ring. 

In summary, we have studied the influence of period- 
ically alternating magnetic fields on the dynamics of a 
quantum state on a tight-binding ring. Our analytical 
and numerical calculations indicate that the evolution of 
a quantum state can be suppressed through the time- 
periodic flux. The efficiency of the scheme depends on 
not only the system parameters 4>x, <fio and u>, but also 
the state parameters fco and a. Based on the features of 
the dynamics, one can choose an appropriate system to 
freeze a given state with an expected average fidelity. It 
can be also exploited to select and hold a specific wave 
packet among the many-body particles, thus providing 
a means to perform the quantum state engineering. We 
expect that such an observation has applications for in- 
formation processing and quantum device physics. 



the similar result but with different optimal </>a ■ The cor- 
responding numerical result for the threshold frequency 
v c (a) as a function of the width of the wave packet is 
plotted in Fig. M 

Based on the numerical results presented in the two 
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